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ABSTRACT 

[-]Methed A method for determining the variations of the resistivity index (RI) of a 
family of rock samples of complex pore structure with at least two pore networks, as a 
function of the water saturation (Sw), and in the presence of a non-oil conducting fluid. 

[-]The volume fraction (f b f 2 , f 3 ) occupied by each pore network and the distribution of 
the pore throats in the various pore networks are determined for each sample of the 
family. The values of coefficients relating the variation of its electrical resistivity as a 
function of its water saturation (Sw) are experimentally determined on a sample of the 
family used as a reference sample. The resistivity index (RI) of all the samples of the 
family is then determined on the basis of the variation of parameters describing the 
layout of the pore network, and using the values of saki-the coefficients measured on the 
reference sample. The method allows to measure measurement of a continuous 
resistivity index curve which is not obtained, as in conventional techniques, from only a 
limited number of points at capillary equilibrium. 

[-] Applications : petrology for example. 
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BACKGROUND OF THE INVENTION 

FIELD OF THE INVENTION 
Field of the Invention 

The present invention relates to a method for determining the resistivity index RI, as 
a function of the water saturation, of certain rocks of complex porosity. 

The method according to the invention finds applications notably in the spher e of 
hydrocarbon reservoir development. 

Descripti on of the Prior Art 

Evaluation of carbonate reservoirs is a particularly difficult task for petrophysicists, 
who still lack precise knowledge concerning the carrying properties within these porous 
media. In relation to siliclastic rocks, carbonates may be simpler on the mineralogical 
plane, but they are incomparably more complex in terms of pore structure and surface 
properties. The largely biological origin of the sediments, combined with various 
diagenetic processes, leads to complex pore structures that may be very different from 
one reservoir to the next. For many carbonate systems, resistivity data calibrations 
carried out in the laboratory contradict the observations made in the field (anhydrous 
production, capillary pressure) and the direct water saturation measurements performed 
on preserved cores. 
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BACKGROUND OF THE INVENTION 

The prior art is notably defined by the following publications : 

- Bouvier L. et al., Reconciliation of Log and Laboratory Derived Irreducible Water 
5 Saturation in a Double Porosity Reservoir, Advances in Core Evaluation, edited by 

Worthington and Longeron, Gordon and Breach Science Publishers, 

- Dixon J.R. et al. (1990), The Effect of Bimodal Pore Size Distribution on Electrical 
Properties of some Middle Eastern Limestone, Soc. Petr. Eng. 20601, 7 th SPE Middle 
Eastern Oil Show, Bahrain, pp.743-750, 

10 - Fleury M. (1998), £«] "FRIM : a Fast Resistivity Index Measurement Method" £»1, 
Proceedings of the International Symposium of the Society of the Core Analysts, Den 
Hague, 

- Fleury M. et al. (2000) [«1 "Frequency Effect on Resistivity Index Curves Using a 
New Method" £»1, Proceedings of the 41 st Annual SPWLA Symposium, Dallas, 

15 - Moore C.H. (2001) £«] ^Carbonate Reservoirs, Porosity Evaluation and Diagenesis 
in a Sequence Stratigraphic Framework" £»1, Developments in Sedimentology 55, 
Elsevier Editions, 

- Petricola MJ.C. et al. (1995) £«1 "Effect of Microporosity in Carbonates: 
Introduction of a Versatile Saturation Equation^ £»1, Soc. Petr. Eng. 29841, SPE Middle 

20 Eastern Oil Show, Bahrain, pp.607-615, 

- Sen P.N. et al. (1997) £«] ^Resistivity of Partially Saturated Rocks with 
Microporosity" £»], Geophysics, Vol.62, No.2, pp.415-425. 
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Understanding and prediction of the effect of the structure of the pore network, of 
the wettability and of the electrical properties of carbone rocks is a real scientific 
challenge, theoretically as well as experimentally. In fact, correct evaluation of these 
parameters has a major impact on estimation of the oil in place, notably for Middle 
5 Eastern giant fields, because the difference in relation to the standard values of Archie's 
exponents m and n in relation to value 2 is so great that estimation of the water 
saturation can vary by more than 20 %. 

Various experimental observations have shown the existence of distinct pore 
populations : micropores, macropores and mesopores, with different coexistence 

10 degrees. In general, the resistivity index curve RI(Sw)=Rt(Sw)/Ro, where Rt is the 
resistivity of the rock to a water saturation Sw, and Ro the resistivity for Sw=l, cannot 
be described by a power law (second Archie's law RI=Sw n ), ire^ that is n is a function 
of the saturation itself. Microporosity can act as a parallel path for the current, which 
leads to a decrease in the values of n and, therefore, to a gradual insensitivity of the 

15 resistivity to saturation, as observed on clayey sands. Microporosity can also be the 
cause of the low measured values of n (typically 1.45). It has also been observed that n 
can increase considerably under certain conditions, and there might be a connection 
between curve RI(Sw) and the capillary pressure curve. The increase in the values of n 
is also a known effect of the wettability, which tends to feveur- favor aqueous phase 

20 discontinuity and therefore to increase the resistivity to water wettability. The effect of 
the wettability can lead to either a sudden increase of n, or to a high value n without 
discontinuity, which can lead to confusion. 



The experimental curves RI(Sw) which have been drawn from various prior works, 
already mentioned above, and from our own observations can have (Fig.l) four distinct 
shapes which do not always meet Archie's laws : 

type I : can be typical of carbonates from the Thamama formation, 

type II : straightens at intermediate saturation and flattens at low saturation (present 
study), 

- type III : single slope at low saturation, extrapolation at Sw=l above Ir=l, and 

type IV : typical of oil-wet systems, high values of n that can increase further at low 
saturation. This is also valid for clasts. 

It can thus be seen that, in the presence of rocks of complex porosity which do not 
meet Archie's laws, a large number of costly measurements is necessary to take account 
of the variability of the pore structure. 

SUMMARY OF THE INVENTION 

The obj e ct of the method according to the invention i&4e-determines the variations 
of the resistivity index (RI) of a family of rock samples of complex pore structure as a 
function of the water saturation (Sw), in the presence of a non-conducting fluid. fr- The 
method of the invention comprises the following stages : 

for each sample of the family comprising at least a first and a second pore network, 
determining the-a_volume fraction occupied by each pore network by applying to the 
various samples an NMR type relaxometry technique!,!; 

- for each sample of the family, measuring by mercury injection the pore throat 
distribution in the various pore networks!,]; 



* 
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- determining experimentally on a sample at least of the family used as thereference 
the values of coefficients relating the variation of ks- the sampled electrical resistivity as 
a function of rts -the sample's water saturation!,]; and 

- determining the resistivity index (RI) of all the samples of the family on the basis of 
5 the variation of parameters describing the layout of the pore network and using the 

values of the coefficients measured on the reference sample. 

According to a first implementation mode suited for samples comprising two pore 
networks, the values of the coefficients relating the-a_total conductivity of the sample to 
the~a_conductivity of the first and second pore networks and to tbe-respective water 
10 saturations of the two networks are determined from the reference sample]^ _the-The 
resistivity index betftg-is_calculated from the respective volume fractions of the two pore 
networks and from the value of the mean saturation from which the network with the 
smaller pores is invaded by the non-conducting fluid. 

According to a second implementation mode suited for samples comprising a third 
15 pore network, the values of the-coefficients relating the»a_total conductivity of tbe-a 
sample to the conductivity of the first and second pore networks and to the respective 
water saturations of the first two networks are determined from the reference sample]^ 
the- The resistivity index bera^is_calculated from the respective volume fractions of the 
three pore networks, from the-a_value of the-amean saturation from which the network 
20 having the smaller pores among the first two pore networks is invaded by the non- 
conducting fluid, and from the value of the mean saturation from which the network 
having the larger pores among the first two pore networks is invaded by the non- 
conducting fluid. 



The resistivity index variations obtained by applying the method are in best 
agreement with the experimental curves whose examples are illustrated in Fig.l. The 
method is based on a powerful experimental technique allowing te-measure ment of a 
continuous resistivity index curve, unlike the conventional technique which describes 
5 curve RI only by a limited number of points at capillary equilibrium. 

BRIEF DESCRIPTION OF THE FIGURES DRA WINGS 

Other features and advantages of the method and of the device according to the 
invention will be clear from reading the description hereafter, with reference to the 
accompanying drawings wherein : 

10 - Figure Fig. 1 shows a diagrammatic representation, in logarithmic coordinates, of 
various types of curves RI as a function of the water saturation Sw observed 
experimentally [;]with curves I to III are -being due to the-pore structure^ ] and curve IV 
is typical of the effects of the-wettability and4 ^which is not specific to carbonates!,]; 

- Figure Fig. 2 shows a pore size distribution of three carbonates studied : (reservoir 
15 carbonate, Estaillade limestone and Brauvillier limestone) obtained by NMR[ ;l~4he 

figur e s indicate with the pore volume fraction being to the left of the vertical 
segmental; 

Figures 3a, 3b Figs. 3a and 3b show the sensitivity of the model with two pore 
networks (DPC) to the saturation Sc at which the micropores are invaded by oil (3a) and 
20 to oc M which is the initial conductivity ratio of the two pore populations (3b)£J; 

- Figures 4 a, 4 b Figs. 4a and 4b show the application of model DPC to sample RC at 
ambient temperature (4a) and under the reservoir conditions with dead oil (4b) £;]w|th 
the experimental curve is- being drawn in- as a thin line and the model m- being drawn as 



a_thick line [; ]with parameters n,, n 2 and a are -being adjusted, f, and Sc afe -being 
measured£,k 

- Figur e Fig. 5 shows an invasion scenario in the model with three pore networks 
(TPC)Ll; 

5 - Figur e Fig. 6 shows an electrical layout modelling of the three networks in the TPC 
model£,]i 

- Figur e Fig. 7 shows a measured (thin line) and model (thick line) curve RI for 
sample EL (Estaillade limestone)[,k 

Figure Fig. 8 shows a measured (thin line) and model (thick line) curve RI for 
10 sample BL (Brauvillier limestone)[J; and 

F igu r e Fig. 9 shows a simplified calibration sequence for the resistivity logs [; ]with 
for non- Archie type rocks, Sw(RI) is necessary [; ] and the laboratory data generally 
provide providing relation RI(Sw). 

DETAILED DESCRIPTION 

15 DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 

W e hav e a c e rtain A number of samples of the same family and we want to 

calculat e exist for which calculation of the resistivity index of each-one sample is 

desired . The method described hereafter allows these indices to be obtained, notably 

from : 

20 - experimental measurements performed on each sample of the family, giving the 
number of pore networks of different size and the volume fraction occupied by each 
pore network, as well as the pore throat distribution for each one, and 
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- an experimental measurement performed on a single sample of the family, used as a 
reference sample, of the variation of its electrical resistivity as a function of its water 
saturation (Sw). 

Nomenclature 
5 Cto = total conductivity at Sw=l 

Cj = conductivity of the i-th population at Swi=l 

fi = NMR volume fraction of pores of the i-th population, i=l,2,3 

RI = resistivity index, Rt/RoSw = mean water saturation 

SWj = water saturation of the i-th pore population 
10 Sc = mean saturation at which the micropores are invaded 

Sm = mean saturation at which the macropores are invaded 

a u = conductivity ratio C x /C 2 and Ci/C 3 respectively. 

Experimental data on samples 

NMR relaxation time measurements 

15 The various pore networks within the rock samples and the volume fraction 

occupied by each network are first determined by means of an -a well-known NMR type 
relaxometry method well known to specialists . An example of such an NMR 
relaxometry method is for example described in French patent application EN-02/... 
filed by the applicant. 

20 We present the measurements performed on three selected carbonate samples. By 

applying an NMR type known relaxometry method known to sp e cialists , w e can see it is 
seen that all the samples are characterized by a double or triple pore size distribution. 
The first one is a reservoir carbonate sample from the Middle-East, classified as a 



packstone. It was cleaned by means of different solvents at high temperature before 
being subjected to the measurements. Under ambient conditions, the experiment was 
carried out with refined oil and synthetic reservoir brine. The experiment carried out 
under reservoir conditions was conducted with filtered degassed reservoir oil at the 

5 reservoir pressure. After cleaning, the wettability of the sample is characterized by a 
moderate water wettability and, after aging, by a high oil wettability. The other two 
samples are outcrop carbonates (water-wet). The porosity of the Brauvillier limestone 
(BL) is essentially intergranular and due to the oolite cortex. The porosity of the 
Estaillade limestone is both inter and intragranular. The pore populations are separated 

10 by a factor 10 at least. 

Resistivity index measurements 

The resistivity index of the samples is then measured. The fast resistivity index 
measurement method (FRIM) described in patent EP-974,839 filed by the applicant 
assignee can be advantageously used. 

15 According to this method, a forced oil-water displacement is carried out on a small 

core of length 2.5 cm and diameter 4 cm for example. This displacement is very close to 
a porous plate displacement method, except that capillary equilibrium is not necessary. 
Under ambient conditions, only two or three pressure stages are used. Under reservoir 
conditions, when dead or crude oil is used, each pressure stage is prolonged for the time 

20 required to obtain wettability stabilization at a given saturation. Drainage is thus carried 
out in approximately two to three weeks, which typically corresponds to the kinetics of 
the chemical processes linked with aging. The simplicity of this method lies in the fact 
that it just requires real-time recording and plotting of the mean saturation and 



10 



resistance to obtain a continuous resistivity index curve free of any artefact. The key 
point is due to the fact that the radial geometry of the electrodes allows to examine 
examining the whole volume of the sample and to compensate for the non-uniform 
saturation profile which appears in the absence of capillary equilibrium. 

5 Under ambient conditions, we us e d the radial-electrode cell described for example 

in the aforementioned patent EP-A-974,839 is used . The complex impedance 
measurements were carried out at a fixed frequency of 1 kHz and the real part was 
extracted to calculate the resistivity index. The highest capillary pressure applied was 12 
bars (for an interfacial tension y=35 mN/m ; under reservoir conditions, with crude oil, 

10 the maximum capillary pressure was reduced proportionally to the interfacial tension). 

Conductivity models 

Double-porosity conductivity model (DPC) 

The goal is to explain the bending of curves RI at low saturation (type I, Figure 
Fig. 1). The double-porosity model considered here is basically very close to those 
15 proposed for clayey sands, where the clays present at the surface of the pores constitute 
a parallel path for the current. For carbonates, we assume the existence of two pore 
networks having parallel electrical conductivities is assumed . The two main ingredients 
in euf ^the model are the description of the pore network invasion during drainage and 
the description of the electrical layout of the two different pore populations. 

20 The first network 1 (macropores for example) represents the major part of the pore 

volume, whereas the second network 2 (micropores) only represents a small fraction 
thereof, which is not necessarily greater than the percolation thresholds. W e first 
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consider th e The saturation of each network, Sv/ X and Sw 2 , which are related to the 
measured mean saturation Sw is first considered by : 

Sw = f^Sw, + f 2 Sw 2 where f, + f 2 =l (1). 

The quantities £ and f 2 represent the pore volume fraction of each population. These 
5 fractions were evaluated by means of the NMR relaxometry technique. We then assume 
that the networks are invaded by oil at different capillary pressures ; the oil reaches the 
small pores at a higher pressure than the pressure observed for the larger pores. From 
the capillary pressure curve previously drawn, this pressure corresponds to a mean 
saturation Sc which can be deduced from the curves obtained by mercury injection 
10 (showing the pore throat distribution). Sw, can be expressed as a function of Sw at high 
saturation Sw : 



Sw, = SW + /l \ sw 2 =1 for Sw>Sc (2) 

Below Sc, relations Sw,=f(Sw) and Sw 2 =f(Sw) require other hypotheses. We-The 
assumption a ss*mre-(0 is_a linear relation and (ii) that Sw, -»0, Sw 2 ->0 when Sw -»0. 

15 Wc deduc e It is deduced that : 
Sw, — Sw, 

ASc 

Sw 2 =— for Sw<Sc (3) 
Sc 
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We now consider th e The conductivity of each network is now considered . In case 
of an initial brine saturation, the total conductivity Cto for the two networks in parallel 
will be : 

Cto = C, + C 2 = Ci (1 + a) where C 2 = a d (4). 

Parameter a is the conductivity ratio of the two 100 % saturated networks. From 
Archie's first law, w e assume that d and C 2 are assumed to be related to the pore 
volume fraction of each population and, therefore, that a is of the order of ratio f 2 /f,. 
When the two networks are invaded by oil, we assum e that each conductivity is 
assumed to be related to the saturation by a power law (as in Archie's second law). The 
total conductivity in the two saturation ranges is as follows : 

Ct = Swj nl d + C 2 forSw>Sc (5) 

Ct = Sw! nl Ci + Sw 2 n2 C 2 for Sw < Sc (6). 



RI = ^ = Sw7 nl 1 + a : for Sw > Sc (7) 

Ct 1 1 + aSV 



RI = Sw; n ' — — ; o _ n2 for Sw < Sc (8) 



1-f-q 
\ + aSw; nX ISwl n2 

By means of equations 4, 5 and 6, we find the resistivity index RI is_as follows : 



When a pore population is in the dominant position (network 1 in the present case), 
functions RI are governed by the resistivity properties (nj) of this population. Equations 
7 and 8 are essentially similar to the formulas used for clayey sands, except that we 
introduce a second exponent n 2 isintroduced which characterizes the second network. 



In the DPC model, there are in total 4 parameters n,, n 2 , a and Sc. For a given 
experimental curve, Sc is measured separately during a mercury injection experiment, 
whereas the other parameters are adjusted. A variation range is however observed for a 
around f 2 /f,, for which a physical explanation exists. At high saturation Sw, the slope on 
5 the log-log scale of RI(Sw) is -n, and, at low saturation, the slope is -n 2 . To a certain 
extent, Sc and a compensate for one another ( Figur e Fig. 3), but a is the most sensitive 
parameter which controls the final value of RI. On e It can be n oted that the 
case Sc=0.05 presented in Figure Fig. 3 corresponds to a situation where the second 
network is not invaded by oil, which gives a horizontal asymptote, 

10 The parameters of the model were adjusted in order to correspond to the 

experimental curve RI measured on sample RC. The value of Sc=0.4 was deduced from 
the elbow of the capillary pressure curve obtained by mercury injection and the volume 
fraction £=0.88 macropores (network 1) was deduced by NMR ( Fig ur e Fig. 2). Despite 
the low permeability of the sample, the lowest saturation obtained is very low (4 %), 

15 which allows good determination of the parameters of the model under ambient 
conditions. The initial slope of curve RI is n^l.71 ; it characterizes network 1. 
Network 2 is very weakly sensitive to saturation (n 2 =0.25) and does not behave like a 
standard network. However, the initial conductivity ratio cc=0.054 of network 1 to 2 is 
of the order of f 2 /fi=0.136. Under reservoir conditions (oil wettability conditions), a high 

20 increase of n x is observed, but curve RI is still not linear on the log-log scale. 
Characteristics n 2 =0.78 and a=0.11 of the second network are also slightly modified, 
but the precision concerning these parameters is lower that under ambient conditions 
because the saturation reached (at the same capillary pressure) is much higher and close 
to Sc. 
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Approximate formula in the case of the DPC model 

Wo propose an An approximate formula reducing relations 7 and 8 to a single 
formula, which is valid for the type I curves is as follows . Wo first note Note that : 

Sw 2 n2 Swf nl « Sc n2 Sw n2 " nl * Sw n2 " nl 

5 with Sw, = Sw and Sc n2 « 1 (15) 

which gives f 2 «fi- Since n 2 «n,, equations 7 and 8 can be approximated by : 

RI = SW nl 1 + C w2 nX (16). 

As can b e sc o n in Figure 10, using Usmgjhe 3 parameters n l5 n 2 and C is sufficient 
to describe the data with precision by means of equation 16. Here, the meaning of C is 
approximately the same as a and it can depend on the temperature. The formulation 
10 given in equation 16 can be used to describe standard experiments carried out at 
capillary equilibrium where only a limited number of points is available. 

Triple-porosity conductivity model (TPC) 

The goal is to explain the straightening of curves RI at intermediate or high 
saturation, and their bending at low saturation on the log-log scale (types II and III, 
15 Figure Fig. l). Wo use here thr > general idea according to which compl e x Complex 
carbonates can have three pore populations referred to, for simplification reasons, as 
micro, macro and mesopores (3, 2 and 1 respectively). As in the double-porosity model, 
we consider the saturation of the 3 populations is considered : 

Sw = fiSwi + f 2 Sw 2 + f 3 Sw3 where f i + f 2 + f 3 = 1 (9) 
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Here again, wo assume that the invasion of these populations by oil during drainage 
is assumed to be sequential. If network 1 is invaded first, we d e fine a mean 
saturation Sm is defined at which network 2 is invaded : 



A 

Sw 2 =1 Sw 3 = 1 for Sw > Sm (10) 

5 Below Sm, it is possible to imagine many scenarios. In general, w e assume linear 

relations for functions Sw,(Sw) and Sw->(Sw) are assumed . According to a possible 
scenario, we assum e that Sw, ->0£ J and Sw 2 ->0 when Sw-^Sc is assumed . Sc is the 
saturation at which the micropores are invaded. It follows therefrom that : 



f x Sm - Sc 

Sw 3 =1 for Sc<Sw< Sm (11) 



10 This invasion scenario is summarized in Figure Fig. 5. Typically, one thinks of a 

situation is considered where f 3 «fi and f,«f 2 , and the micropores (network 3) are 
invaded at a pressure that is much too high to be observed during the experiment. 

Consider now the conductivity of these populations. W o assum e that n e tworks 
Networks 1 and 2 are assumed in series, whereas network 3 is in parallel as shown in 
15 Figure 6. The arrangement in series may seem to contradict a sequential invasion by oil 
(which is an essentially parallel mechanism), but such a situation is possible. At Sw=l, 
conductivity Cto of the system shown in Figur e Fig. 6 is as follows : 

Cto = (Cf 1 + C2" 1 )" 1 + C 3 = Q [(1 + <xf V + ot 2 ] 

where C 2 = oci Ci, C 3 = ot 2 Cj (12). 
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In the two saturations ranges, wo obtain i t is obtained that : 

Ct = Ci [(Swf nl + af 1 )" 1 + a 2 ] for Sw > Sm (13) 
Ct = Ci [Swf nl + Sw 2 " n2 af 1 )" 1 + a 2 ] for Sc < Sw < Sm (14). 

The resistivity index can be calculated from equations 12 to 14. For the measured 
5 values off,, f 2 and f 3 , n,, n 2 , a! and a 2 have to be adjusted to the experimental data. As in 
the DPC model, a! is expected to be of the order of f 2 /f, and cc 2 of the order of f 3 /fi. 

This model was tested with resistivity index curves measured on samples EL and 
BL. On sample EL (Figure 7), straightening of the curve at high saturation is 
reproduced qualitatively at high saturation and flattening at low saturation is 

10 reproduced. Wo observe a A discontinuity at Sm (deduced from the elbow of the 
capillary pressure curve obtained by mercury injection) is observed b ecause of a sudden 
modification of relation Sw,(Sw) (as shown in Figure Fig. 5). After adjustment, the two 
dominant pore populations have similar characteristics (n,=n 2 =1.5). The third 
population, which behaves like a parallel circuit, is not invaded by oil and it is only 

15 characterized by conductivity ratio a 2 . On sample BL ( Figure Fig. 8), curve RI can be 
explained by a similar mechanism of resistance in series, but with a very high value Sm. 
Wo finally observ e that th e Thetwo networks observed have very different values n and 
that network 2 is much less conducting than network 1 (cti=7.1). 

The method according to the invention finds applications notably in the sphere of 
20 hydrocarbon reservoir development. 



